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Motivated by recent experiments, we consider theoretically the compression of droplets pinned at 
the bottom on a surface of finite area. We show that if the droplet is sufBcicntly compressed at the 
top by a surface, it will always develop a shape instability at a critical compression. When the top 
surface is flat, the shape instability occurs precisely when the apparent contact angle of the droplet at 
the pinned surface is tt, regardless of the contact angle of the upper surface, reminiscent of past work 
on liquid bridges and sessile droplets as first observed by Plateau. After the critical compression, 
the droplet transitions from a symmetric to an asymmetric shape. The force required to deform the 
droplet peaks at the critical point then progressively decreases indicative of catastrophic buckling. 
We characterize the transition in droplet shape using illustrative examples in two dimensions followed 
by perturbative analysis as well as numerical simulation in three dimensions. When the upper surface 
is not flat, the simple apparent contact angle criterion no longer holds, and a detailed stability 
analysis is carried out to predict the critical compression. 



I. INTRODUCTION 

The interaction between liquids and solids is ubiqui- 
tous in our daily life, from droplets on a windshield to 
ink in our printers [1]. Capillary phenomena arise as a 
consequence of intermolecular forces and manifest them- 
selves on large scales by the tendency of liquids to min- 
imize their surface area [2]. The theory for the shapes 
of droplets was proposed over two hundred years ago by 
Young [3] and Laplace [4], and since that time much has 
been learned on the wetting of solids by liquids [5, 6]. Of 
interest to us in this paper is why certain droplet config- 
urations may be unstable. Plateau observed that a liquid 
jet would ultimately break up into droplets because the 
energy of the initial cylindrical shape is lowered by long 
wavelength perturbations [7]. The dynamics of this in- 
stability was later elucidated by Lord Rayleigh [8] while 
more recent work has explored the finite-time singularity 
at the break-up [9]. 

A somewhat lesser known observation made by 
Plateau, in the same work [7], details how a droplet 
suspended between two equal circular disks (a so-called 
liquid bridge) loses axisymmetry when sufficiently com- 
pressed. Plateau [7] observed that the onset of this shape 
instability occurs when the profile of the droplet at the 
point of contact with the pinned surface becomes tan- 
gent to the disks, hence making for an apparent droplet 
contact angle of tt (see also the review in Ref. [10]). An 
analytical solution for the shape of a liquid bridge was 
put forward by Howe [11] along with an initial analy- 
sis of the instability [11, 12]. It was later shown formally 
that if the axisymmetric droplet shape is described in po- 
lar coordinates by a single valued radius function then it 
is always stable to asymmetric perturbations for pinned 
boundary conditions [13]. Such single- valued description 
breaks down when the droplet becomes tangent to the 
disks. Russo and Steen showed then that past this point 
the droplet is unstable to asymmetric perturbations [14] . 

Subsequent works [15-18] further elucidated the space 



of stability of these liquid bridges for which the results 
of Plateau is only one of many. In particular, it was sim- 
ilarly shown that a pinned droplet deformed by a gravi- 
tational field, rather than by compression from an upper 
surface, will also transition to an asymmetric shape past 
the point when the droplet profile is tangent to the pin- 
ning line at the point of contact [10]. A related shape 
instability arises when two free-surfaces are squeezed to- 
gether [19]. If two droplets, or bubbles [20, 21], are 
brought into contact, then initially the interface sepa- 
rating them is perpendicular to the direction of compres- 
sion and their shapes deform in an axisymmetric fashion. 
Past a critical conformation the separating interface ro- 
tates and the droplets, or bubbles, lose axisymmetry. A 
general framework for the stability of equilibrium states 
of capillary phenomena is provided in the text by Myshkis 
et al. [22]. 

In recent experiments by Nagy and Neitzel, droplets 
pinned to a bottom surface and compressed by a per- 
fectly non-wetting surface from above have been shown 
to develop a geometric asymmetry at a critical defor- 
mation [23]. Experimentally, the perfectly non- wetting 
condition on the top surface is obtained using a thin 
layer of air maintained by thermocapillary convection be- 
tween the cold surface (of arbitrary contact angle) and 
the hot droplet, leading to an effective contact angle of 
180° [24, 25]. On the bottom surface, pinned boundary 
conditions are achieved by extruding liquid through a 
capillary, but one can also imagine a droplet confined to 
the top of a small post or a disk. When the droplet 
is compressed, by displacing the top surface, below a 
threshold the droplet maintains axisymmetry, but upon 
reaching a critical conformation, the droplet bulges to 
one side indicating a shape instability. Subsequent Sur- 
face Evolver [26] simulations of the same setup verified 
the geometric asymmetry [23]. 

In this paper we show that the shape instability ob- 
served by Nagy and Neitzel [23] arises with the same ge- 
ometrical criterion as that of the Plateau liquid bridge. 
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FIG. 1. Images of a water droplet at the tip of a straw compressed against a liydrophobic surface; (a): The droplet after initial 
contact, with pinned boundary conditions at the bottom and a fixed hydrophobic contact angle at the top; (b): Axisymmetric 
deformation of the droplet prior to the onset of shape instability; (c) Asymmetric bulge-like shape of the droplet after the 
critical compression. 



namely when the droplet surface becomes tangent to the 
pinned surface at the point of contact (apparent contact 
angle of tt). We demonstrate that this instability crite- 
rion does not depend on the droplet contact angle on the 
upper surface. We reveal, however, that that the geomet- 
ric instability criterion no longer holds when the deform- 
ing surface is not flat, as is illustrated with a compression 
of the droplet by a conical surface. We first use a simple 
table-top experimental example to show that the shape 
instability occurs even if the deforming surface is not per- 
fectly non-wetting. We then employ a two-dimensional 
analysis to provide a preliminary analytical approach to 
the instability, and derive the stability criterion and its 
independence on the top contact angle. We next carry 
out a three-dimensional perturbation energetic analysis 
confirming the two-dimensional theory, and demonstrat- 
ing the breakdown of the geometric stability criterion 
when the deforming surface is no longer flat. We finally 
utilize numerical computations using Surface Evolver to 
confirm our theoretical predictions. 

II. TABLE-TOP EXPERIMENTS 

We first performed simple table-top experiments in an 
effort to observe whether the shape instability occurs in 
the case where the compressing surface is not perfectly 
non- wetting. The results are shown in Fig. 1. We use a 
plastic straw extruding water, coated with hydrophobic 
spray to delay its wetting prior to the shape instabil- 
ity. The water droplet is put in contact with a slightly 
hydrophobic surface. In Fig. 1 we see that initially the 
droplet is symmetric (a), and remains symmetric under 
small deformations (b). When the droplet has been suf- 
ficiently squeezed, it rapidly develops a pronounced geo- 
metric asymmetry (c). 

III. ENERGETIC ANALYSIS 

Consider a droplet pinned along one surface, which 
we will call the bottom surface, along a constant circu- 



lar area. The droplet has an volume V which remains 
constant by conservation of mass. We assume that the 
droplet is smaller than the capillary length and hence 
ignore gravity [1]. The droplet shape is the one which 
minimizes the functional 

[ jdS- [ pdV, (1) 
Js Jv 

where 7 is the interfacial tension coefficient which is in- 
tegrated over all interfaces S and the pressure p is a La- 
grange multiplier enforcing mass conservation. We fur- 
ther assume that the droplet is in contact by a second 
surface, parallel to the bottom surface, and which we 
will refer to as the top surface. That surface, of area W, 
does not pin the drop but instead the droplet contact 
line is free to move along the top surface with the set 
contact angle 9. We ignore effects arising from contact 
angle hysteresis. In that case the relevant surface energy 
now becomes 

E = "fsv{W - Ssi) + jsiSsi + Siv^iy + const., (2) 

where the subscripts sv, si and Iv indicate the solid- 
vapor, solid-liquid and liquid-vapor interfaces respec- 
tively. The surface energy of the pinned interface is un- 
changing and thus rolled into the constant term. Using 
Young's equation for the equilibrium contact angle 9 on 
the top surface [1], Eq. (2) may then be recast as 

E = jiv (Sly - Ssi cos 9) + ^svW + const. (3) 

Given that W is constant, we see that finding an equi- 
librium is tantamount to minimizing the projected sur- 
face S = Sly — Ssi COS 9. Note that when the solid sur- 
face is nonwetting, 9 = tt, and one has thus to minimize 
S — Ssi + Sly. From this point on we neglect the con- 
stant terms in Eq. (3), and drop the Iv subscript for the 
liquid- vapor interfacial coefficient for simplicity, hence we 
denote E = jS. 
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FIG. 2. Schematic representation of the droplet in various stages of deformation, for a contact angle — 27r/3. (a): No 
deformation; (b): Droplet at the critical conformation, 7? — R^, (c): Asymmetric droplet post threshold R = Ra {S > 0); (d); 
Symmetric droplet post threshold R = Rs {5 > 0). The droplet volume (termed A in 2D, V in 3D) remains constant; the initial 
radius Ri and surface area Si are thus set by the value of the droplet volume. 



A. Two-dimensional analysis 

1. Flat surface 

In order to get a sense of the essential physics of the 
droplet instability we first illustrate the process in two 
dimensions. The geometric criterion for the instability 
of symmetric solutions is easier to derive and visualize 
than in three-dimensions, which is a calculation we will 
address in the next section. 

The setup and notations are illustrated in Fig. 2. The 
volume of the droplet is denoted V in three dimensions 
and A in two dimensions, and the interfaces are denoted 
S in both cases. The liquid-vapor interface which mini- 
mizes the total interfacial energy is a surface of constant 
mean curvature, and thus in 2D it is circular. We there- 
fore write the shape of the liquid-vapor interface, param- 
eterized by t in the {r,z} plane, as 

r{t) = ro + Rcos{t), 

z{t) ^ zo + Rsm{t), (4) 

where R is the radius of curvature of the 2D droplet 
originating from the center {rQ,zo), and t is the polar 
angle about that point measured positive counterclock- 
wise from the point where the droplet is locally parallel 
to the z axis on the right-hand side. The contact point 
where the droplet is pinned occurs at t = ti, and the 
point of contact with the wall W occurs at t = i2- We 
define At = t2 — ti so that the surface area of the liquid- 
vapor portion of the droplet is Siy = 2RAt. The contact 
area with the bottom surface is assumed to be of length 
L = 2ri . We subsequently employ the subscript 1 or 2 on 



any function of t to denote the function evaluated at ti or 
t2 respectively. With this in mind, the height between the 
two surfaces (clearance) is given hy h = Z2 — zi. We use 9 
to denote the equilibrium contact angle of the droplet on 
the top surface, and /3 the apparent contact angle with 
the bottom, horizontal, surface (where the contact line is 
pinned) . 

When the top surface is not in contact with the 
droplet, as in Fig. 2a, then t2 is taken to lie on the 
axis of symmetry of the undeformed droplet {t2 — 7r/2). 
The initial droplet volume, which remains constant, is 
A = (i?2/2)(2At + sin[27r-2At]), and the pinning length 
i2 = 4i?f(l - sin^[(7r - 2At)/2]) closes the system. The 
droplet will form a half circle when L = 4^/ A/2Tr. 

When there is contact with the wall, the droplet takes 
the shape illustrated in Fig. 2b, and two 'disconnected' 
liquid- vapor surfaces are present. In order to maintain 
mechanical equilibrium these two surfaces must have an 
identical radius of curvature, R. The point t2 is located 
where the droplet is in contact with the surface, and 
therefore we necessarily have t2 = 6* — 7r/2. If the appar- 
ent contact angle at the pinned line, /3, satisfies /3 < tt, 
in other words ti > —it/2, we will show below that the 
only solution for the two-dimensional droplet is left-right 
symmetric (the details below will demonstrate that there 
exists no asymmetric solution). 

At the critical point (denoted with a subscript c), 
which is depicted in Fig. 2b, we have ti = t\ = — 7r/2 
and At = 9. The total surface energy is then given 
by Ec/-/ = {2eRc - Lcos9), where = 9- cos 6* sin 6*. 
Here the radius Rc is given by conservation of volume, 
A = QR^ + LRc{l — COS0); this quadratic equation for 
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i?c can be solved to obtain the critical radius for a given 
initial droplet volume, pinning length and contact angle. 
The critical gap height, /i^ = Z2 — = Rc{l — con 6), 
increases monotonically with the contact angle (which 
ranges from zero to ^ = tt). 

If we now reduce the gap height by some amount 
S, so that h = h'^ — 6, two possibilities exist for the 
droplet shape, either asymmetric (Fig. 2c) or symmet- 
ric (Fig. 2d). In each case the left and right sides of the 
droplet need to have the same radius of curvature to en- 
sure mechanical stability. Prom this point on all lengths 
will be non-dimensionalized by the critical radius Rc, and 
hence i?c = 1- 

In the asymmetric case, denoted by the subscript a, 
one of the circular ends moves in and the other moves 
out (see Fig. 2c). The droplet volume in this instance is 
A = QRl + hL, meaning that Rl = 1 + LS/Q. If we 
take 5 < Q we see that Ra < 1 which is impossible as 
the droplet would not span the required length h, hence 
there can be no asymmetric shape for S < 0. Further, we 
see that for (5 > 0, i?a > 1- The relevant surface energy 
is given by Ea/^ = ^ORa — LcosO in that case. 

Alternatively, the droplet may remain symmetric (sub- 
script s) past the critical point, as shown in Fig. 2d. 
Conservation of volume leads to A = i?^[G — 7r/2 — ii -|- 
cos ti (sin ti—2 cos 9)]+hL, where we set the gap height to 
heh= 1— cos 9— 6, and hence = — sin~^ [cos 9 + h/Rs]- 
Note that when in contact with hydrophilic surfaces, 
9 < n/2, droplets undergo a inversion of the concavity of 
the liquid- vapor interface for larger negative values of S 
where this geometry no longer holds. The relevant sur- 
face energy is now give by £'5/7 = 2(6 — ti — 'k/2)Rs — 
{L — 2Rs cost]) cos 9. 

The symmetric and asymmetric cases are the only two 
possible mechanically stable solutions. The one which 
appears in equilibrium is the one which minimizes sur- 
face energy. Solving for Rg and Eg, can in general only be 
done numerically. However, the difference in the surface 
energies between the symmetric and asymmetric config- 
urations, AS = Eg — Ea, may be formally bounded as 
follows 



AE 
2^ 



e + A - V02 + e(A + r) > o. 



(5) 



where = 9 — sin cos 0, X ~ a ~ sin a cos and F = 
(cos a — cos^) sin a and while 9 G [0,7r] and a = —{tf + 
tt/2) e [0,^] (the superscript s indicates the symmetric 
contact point). We see that evidently G > and A > 
as 9 > sin0. Rearranging we obtain 



A-r > - 



Substituting for back we find 



a — sin a cos a > — 



6' 



(6) 



(7) 



since the left-hand side is non-negative the inequality is 
proved. The equality holds only at the critical point when 



tf = ti = —n 12. in other words when (5 = and Rg = 1. 
We have thus proven that, past the critical point {5 > 0), 
the asymmetric conformation is always the droplet shape 
which minimizes the free energy, independently of the 
value of the contact angle on the top surface. Droplets 
after this critical compression are thus expected to always 
display an asymmetric shape. 

For small values of 5, the value of Rg can be found by 
hand, and provided the contact angle 9 is not too small 
we obtain 

Rs = l- — + (^[l_-COsg]+29)^ ^, ^ ^^^3), 



1 — cos t 



256sin'"(6'/2) 



(8) 



Eg 



= 2Q-Lcos9 + L5 + 



L{1- cos9) + Q 2 /-o^xs 



4sin''(6'/2) 



(9) 



We then get 



AE _ {L[l- cos 9] + 20) 2 
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16esin''(6'/2) 



+ 0{5^) > 0. (10) 



Notably, the force required to deform the droplet is 
proportional to the spatial rate of change of the surface 
area hence, for 5 > we find 



/ 



dE 

'as 



(11) 



This may be recast in terms of the pressure jump f/L = 
Ap = j/Ra which recovers the Young-Laplace equa- 
tion in two dimensions. Since the radius of curvature 
has increased, the pressure has decreased in this con- 
figuration. The slope of the force at the critical point 
df /d5\g^Q = — 7L^/29. As might be physically expected 
the pressure drop is least steep for large L and 9 = tt. 

The two-dimensional results are illustrated in Fig. 3 
where we plot the surface energies, E (main plot), and 
the energy slope (force, inset) against the change in gap 
height from the critical conformation, S (for the values 
L = 2 and 9 = 27r/3). The surface energies increase 
monotonically as the droplet is progressively compressed. 
Beyond the transition point {6 = 0) both the asymmet- 
ric conformation (black) and the symmetric conforma- 
tion (gray) are geometrically permissible although the 
asymmetric conformation is always energetically favor- 
able. Note that the solid gray line is obtained numeri- 
cally for the symmetric conformation whereas the dashed 
line is the asymptotic result. 

As seen above, the slope of the energy is the force re- 
quired to deform the droplet. We therefore see that in- 
creasingly more force is required to compress the droplet 
up to the critical conformation, but that beyond the tran- 
sition point the more we compress the droplet the less 
force is required to deform it. The non-monotonic force 
profile is reminiscent of the non-monotonic pressure re- 
quired to infiate a balloon [27]. For squeezed droplets. 
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FIG. 3. Surface energy, E, vs. a change in separation, 6, witli 
L — 2 and 6 = 2tt/3 (two-dimensional case). The solid lines 
are exact solutions and the dashed line is asymptotic. The 
inset shows the slope of E. We see that when 5 — (indi- 
cated by the dotted line) the surface energy bifurcates into 
symmetric and asymmetric conformations. The asymmetric 
conformation (solid black line) has lower surface energy than 
the asymmetric conformation (solid gray line) and hence is 
the energetically preferable state. 

that result means that the droplet will buckle at the tran- 
sition point if an increasing (or constant) load is applied. 
Unlike the classic Euler-buckling of beams [28] which 
can support an increasing load after a buckling event, 
the buckling instability of a droplet is catastrophic, a so- 
called limit-point instability, as the droplet can no longer 
sustain the force at the point of instability and collapses. 
In the context of superhydrophobic surfaces in which a 
droplet rests on a series of posts this would be what has 
been referred to as the impalement transition between 
the Cassie and Wenzel states [6, 29]. Precisely the same 
behavior will be obtained in three dimensions. 



2. Inclined surface 

How different is the transition to an asymmetric shape 
is the upper surface is not flat? Since we do not want 
to embed a broken-symmetry in the top surface we insist 
that the top surface remains axisymmetric, the simplest 
example of which is a cone of slope x ~ or a wedge in 
two dimensions. One finds that a droplet undergoes a 
similar symmetric to asymmetric bifurcation of possible 
solutions precisely at the moment when the droplet on 
the bottom surface (at t = ti) is parallel to the upper 
surface, in other words, when the apparent contact an- 
gle with the bottom (pinned) surface reaches the value 
P = TT — X- As we show below, this two-dimensional 
criterion is no longer indicative of the onset of the insta- 
bility in three dimensions where the critical compression 
becomes a more complex function of the droplet volume 
and contact angle. The quantitative agreement between 
two and three dimensions is therefore restricted to the 
case where the top surface is flat. 



FIG. 4. Schematic representation of an axisymmetric pinned 
droplet compressed by a cone of slope x- The angle 6 indicates 
the contact angle with the top surface {t — 12) and /3 the 
apparent contact angle with the bottom surface {t = ti). Here 
Vi is the volume of the liquid and Vi is the volume of revolution 
inside the cone from its tip to t2- 

B. Three-dimensional analysis 

1. Axisymmetric solution 

In three dimensions the problem becomes more com- 
plicated due to the introduction of an second radius of 
curvature. Unlike the two-dimensional case, asymmetric 
analytical extremum are unavailable. We follow here the 
approach of Russo and Steen [14] by first assuming that 
the droplet is axisymmetric, which yields analytical so- 
lutions. We then do small asymmetric perturbations to 
the axisymmetric shapes to find the configuration where 
the axisymmetric solution is no longer energy minimiz- 
ing. Finally in Section IV we use Surface Evolver[26] 
simulations to explore the mechanical properties of the 
asymmetric shapes. 

The energy functional, Eq. (1), may be written 

£ = [ jdA- [ ^cosOdA- [ pdV. (12) 

The axisymmetric droplet is schematically depicted in 
Fig. 4. It is easier to integrate V = Vi + Vi and so we have 
to subtract the volume Vi. A known function g{r, z) = 
defines the topology of the top surface. Assuming ax- 
isymmetry and parameterizing the shape by t we may 
recast the integral as 

£ = J TT (^2'-fr\/ + — pr^z^ dt 

-'-/cos0Ssi{r2)+pV,{r2), (13) 

or simply £ = J^^ F{r, f, i, t)dt+f{r2), where again num- 
bered subscripts indicate evaluation at ti or t2 and the 
over-dot represents a partial derivative with respect to 
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t. Note that since the topology of the upper surface is 
known then the only unknown for Ssi and Vi is the loca- 
tion of the contact point, r2. 

Extremizing Eq. (13) leads to the Young-Laplace equa- 
tion 



d 
dt 



2jrz 



■ pr 



0, 



and the boundary condition 

df OF _ dF dg/dr 
dr2 dr dz dg/dz 



0. 



(14) 



(15) 



t=t2 



Eq. (14) was shown by Howe [11], and restated in En- 
glish by Gillette and Dyson [12], to have a solution 



rit) 
z{t) 



= A\/ 1 - sin^ Osin^ t, 
A[E2{n,t) + Ei{n,t) cosf^], 



where 



Jo 



(16a) 
(16b) 



(17) 



are elliptic integrals of the first and second kind. The 
solution described by Eq. (16) has two parameters: the 
first one, A, with units of length, merely acts as a scaling 
factor; the second parameter, D,, which is dimensionless, 
modulates the shape, and may be interpreted through 
the mean curvature H as follows 



cos ft 



\ // 

AH 



(18) 



From this point on wc rescale all distances by A (equiva- 
lently we set A = 1). Note that using these parameters, 
the pressure is then given by 



P 

7 



A(l -Fcosfi) 



(19) 



We note that the energy functional in Eq. (13) may be 
recast using a unit speed parameter s (arclength of the 
generatrix) [22] which then simplifies the formulae shown 
here; however, because the analytical axisymmctric solu- 
tion is nonunit speed[17, 22] we leave the formulation 
general. 



Using Eq. (21) into Eq. (15) we obtain the contact angle 
condition 



Z2 = 7-2 tan(^ + x)- 



(22) 



By differentiating our solution, Eq. (16), and substi- 
tuting into Eq. (22) we can then obtain the location of 
the contact point as 



t2 




arccos 



- cos{9 + x) cot 



(23) 



Note that if we were to select 9 = ty /2 — x then wc ob- 
tain t2 = as expected. Also since the droplet shape is 
symmetric about the point t ~Q then the starting point 
(pinned contact point on the bottom surface) is given by 
the negative of the formula for t2, and we have 




-|- arccos 



■ cos j5 cot 



(24) 



where in Eq. (24) (3 is the apparent contact angle with the 
horizontal bottom surface. The point where /3 = tt is re- 
ferred to as the Steiner limit [14]. In the two-dimensional 
analysis, recall that we obtained that the droplet goes 
unstable when the apparent contact angle, P, is parallel 
to the upper surface and hence ti = tl when /3 = tt — 
We show below that this is true only when the wall is flat 
X = 0, but does not remain valid when x 0. 



3. Physical contact 

If the slope of the cone is nonzero, x Oj we run 
the risk of having physical contact between the upper 
surface and lower surface prior to the instability hence 
we have to limit the regime of f2 and t\ (or /3) to exclude 
this possibility. Since the parameter A merely scales the 
shape it is inconsequential. Solving the equation 



z{t2,^) - z{ti,Vl) = r(t2, ^^) tanx, 



(25) 



yields the limit of physically realizable solutions in the 
Q,-t\ plane. 



2. Boundary conditions 



4- Perturbations 



If the upper surface is a cone with slope x as in Fig. 4 
then it is described by 



g{r,z) = rtanx- {z - h) = 0, 



(20) 



where h is the smallest separation between the pinned 
surface and g. The flat wall is included here as the special 
case X = 0- We now have 



j{r2) = -JcosO hp— tan x- 

cosx 3 



(21) 



In order to observe the instability of axisymmetric 
shapes, following the method of Russo and Steen [14] , we 
add to our axisymmetric solution (dcflned by r{t) and 
z{t)) small, non-axisymmetric perturbations deflned by 
the functions F{t, (j)) and G{t, (p) which contain explicit 
dependence on the azimuthal angle <j). The liquid-vapor 
interface is thus given by the following equations 



i?(i,< 
Z(i,, 



r{t) + F(t,(l)), 
z{t) +G{t,(P). 



(26) 
(27) 
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f2 
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b) 



FIG. 5. The onset of the asymmetric instabiUty with varying shape Q. The data points in this figure indicate the Hmit of 
stability found by solving AE{ti) — for a fiat wall (case a,, x ~ 0) and a conical top surface (case h, x ~ tt/S). The contact 
angle on the top surface varies between 7r/4 and tt by 7r/4 increments. The dashed line (purple online) indicates the point 
when the apparent contact angle (3 — n (droplet tangent to the bottom surface) while the solid line (blue online) indicates 
when P = n — X (bottom contact line parallel to the top surface); these two limits coincide in the fiat wall case of (a). The 
dash-dot lines indicate the limit of physically realizable solutions when X 7^ 0- Our parameter a £ (—0.06, —0.29) is selected 
to give the earliest possible instability. We see when x — (flat wall) the limit of stability corresponds always to /3 = tt, but 
when X 7^ the data does not always lie on either the /3 = 7r — xor/? = 7r curves. The dotted lines (red online) in each plot 
are examples of a paths of constant volume with arrows indicating the direction of decreasing gap separation; they are further 
illustrated in Figs. 7 and 8. 



The perturbations given by F and G must yield shapes 
which satisfy the boundary conditions and preserve vol- 
ume, AVi = Vi{R, Z) — Vi{r, z) = 0, in order to be viable. 

To discover the onset of the instability we solve for 
t'i{Vl) such that surface energy of the perturbed shape 
is equal to the axisymmetric shape /S.E — E{R, Z) — 
E{r, z) = 0. This will yield the point of the instabil- 
ity for any axisymmetric shape fl. The results of this 
perturbative analysis are given below while the mathe- 
matical details are given in Appendix A. The perturba- 
tion analysis can only reveal the onset of the asymmetry; 
to probe the mechanical properties of the asymmetric 
droplets Surface Evolver simulations are used in Section 
IV. Alternatively, other schemes [17, 18], often devel- 
oped for liquid bridges, might be adapted for use with 
the contact angle boundary condition. 

5. Results 

Our theoretical results are illustrated in Fig. 5. We 
plot the dependance of the critical point at which the ax- 
isymmetric state becomes unstable to non-axisymmetric 
perturbations, t1{n), with the shape parameter, f2, for 
contact angles on the upper surface ranging from 7r/4 to 
TT. The case of a flat wall is shown in Fig. 5a while the 
conical surface with x = ""/S is displayed in Fig. 5b. In 
both figures, the dotted lines (red online) are examples 
of a paths of constant volume with arrows indicating the 
direction of decreasing gap separation. 



Plotted in Fig. 5 as a dashed line (purple online) is the 
Steiner limit which corresponds to an apparent contact 
angle on the (bottom) pinned surface of /3 = tt, while 
the solid line (blue online) indicates where the appar- 
ent contact angle at the bottom surface is parallel to the 
top surface, i.e. f3 = ir — x- In the case of a flat wall, 
X = 0, the two lines are coincident. For this geometry 
we see that for any contact angle on the upper surface 
the point of instability occurs on or immediately after 
the Steiner limit. Since we use only a truncated and 
restricted series to represent the asymmetric perturba- 
tions (see Appendix A), this results represents an upper 
bound. Using a symmetrization argument, Gillete and 
Dyson showed that the droplets in a liquid bridge[14- 
18], which can be represented by single valued functions 
r{z) (in other words droplets prior to the Steiner limit), 
are stable to axisymmetric perturbations[13]. There is 
nothing in their argument that prohibits its application 
when we permit r2 to be variable so long as Z2 is fixed; 
however, such a map does not preserve contact angle. 
To apply in our case when x = we must adjust the 
contact angle, but it can easily be shown that one can 
always smoothly change the contact angle with a vanish- 
ing change in the surface area and volume. As lower and 
upper bound coincide, we get the final result that, in the 
fiat wall case, the limit of stability is when the droplet 
is parallel to the pinned surface at the point of contact 
(apparent contact angle of tt), similarly to the stability 
threshold for liquid bridges and sessile droplets [10, 18]. 

The case where the top surface is conical (x — tt/S, 
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Vijrl 

FIG. 6. Critical separation h'^ vs liquid volume Vi 
non-dimensionalized by the pinned radius r\ for Q = 
{■k/4, 7r/2, 37r/4, tt} and x = 0. 

Fig. 5b) shows a different stability behavior. Because 
the upper surface may now contact the lower surface be- 
fore the droplet goes unstable, a portion of parameter 
space, depending on the contact angle, is excluded; this 
is indicated by the dash-dotted lines. Our stability cal- 
culations show, in this case, that the droplet does not 
always go unstable when its surface is parallel to the up- 
per surface (blue solid line) nor when it is tangent to the 
lower siirface (purple dashed line). The two-dimensional 
prediction cannot therefore be extended to three dimen- 
sions. 

Unlike, for the flat surface, the symmetrization argu- 
ment clearly does not hold when x 7^ and hence our 
prediction here for the limit of stability can only be re- 
garded as an upper bound. We will show in the next 
section that this prediction can perform well neverthe- 
less. 



IV. SURFACE EVOLVER COMPUTATIONS 

In order to corroborate our asymptotic predictions we 
appeal to numerical simulations using Surface Evolver 
(SE) [26]. Nagy and Neitzel also used the same program 
to confirm the existence of their observed instability in 
an idealized setting [23]. In Surface Evolver the droplet 
shape is discretized into triangular facets whose positions 
are defined by a position vector X hence the energy is 
i?(X) [30]. At a minimum of energy we have 'VE = 
and the Hessian matrix, H = d'^E/dXidXj, must be pos- 
itive definite. We run SE simulations of both the flat and 
conical geometries and see that, as the separation h ap- 
proaches a critical separation h^, the smallest eigenvalue 
of H approaches zero. If the separation is reduced below 
this critical separation then the symmetric form yields 
negative eigenvalues and the true minimum is then given 
by an asymmetric shape. 

Appealing to dimensional analysis tells us that we 



can write the critical separation in the form h'^/ri = 
^(Vi/rf,9,x) [14, 31]. Given that we have an analytical 
equation for the the limit of stability when the compress- 
ing top wall is flat (x = 0) in Eq. (24), we can plot this 
function easily for various constant 9, with results shown 
in solid lines in Fig. 6. We then run Surface Evolver, the 
symbols in Fig 6 indicate simulations performed where 
h = h'^ was stable but h = 0.99h'^ unstable. This then 
demonstrates that the instability occurs precisely as pre- 
dicted theoretically in the previous section. We note that, 
as predicted by the two-dimensional example, the insta- 
bility is the least physically apparent for large volumes 
and high contact angles, 6 n. In Surface Evolver this 
manifests itself by a need to use very high resolution to 
properly capture the onset of the stability in these phys- 
ical regimes. 

As a more concrete example we define a specific phys- 
ical system and illustrate its transition to asymmetry. 
We choose as an example a droplet of volume V = 4 
pinned on a base of radius ri = 1 (arbitrary units) com- 
pressed by a flat upper surface of contact angle 6 = 120° . 
Solving for T/(A, O, (O)) = 4 and r{A,n,tl{n)) = 1 we 
determine theoretically that A « 1.42 and il w 2.089 at 
a predicted critical separation h'^ « 0.6923. The con- 
stant volume path of this system, from h = l.bh" to h", 
was illustrated in parameter space in Fig. 5 by dots (red 
online), the arrow indicating the direction of increased 
compression of the droplet. 

In Fig. 7 we show the resiilts of Surface Evolver simula- 
tions of this system with the separation, h, ranging l.bh'^ 
to O.Gh". First and foremost, we see that the instabil- 
ity occurs precisely at the analytically predicted h = h'^. 
The pressure in the droplet increases until the critical 
point where the smallest eigenvalue in the system van- 
ishes, indicating a loss of stability. Beyond the critical 
point, symmetric shapes have negative eigenvalues indi- 
cating an unstable saddle. The stable shapes beyond 
the critical point are asymmetric and have progressively 
lower pressure. In Fig. 7b we display the variation of the 
surface energy of the droplet, E, with the separation dis- 
tance between the surface; in inset we plot the slope of 
the energy (obtained by numerical differentiation of E), 
i.e. the force acting on both surfaces resisting compres- 
sion. The surface energy progressively increases as we 
deform the droplet from h = 1.5/i^ to h = 0.6/?/'. The 
slope of the energy (shown in the inset) indicates that 
the force required in order to further deform the droplet 
peaks at the critical separation and then monotonically 
decreases with increasing deformation. 

Experimentally, if we were to progressively load a 
droplet with increasing force then at the point when the 
droplet is tangent to the pinned area it would become 
mechanically unstable, regardless of the contact angle of 
the deforming surface (and true for both two- and three- 
dimensional systems). The system would undergo a dy- 
namic collapse past this point. These mechanical consid- 
erations are important for the design of systems such as 
the one described by Neitzel et al. for load support in a 
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FIG. 7. Surface evolver simulations of the onset of the instabihty for a flat wall, x = with V = 4 and ri — 1. (a): Pressure, 
p (axis on the left, circles joined by a dashed line is numerical data and the solid line is analytical), and smallest eigenvalue 
of the Elessian matrix, Ai (squares with axis on right), versus separation h. We see that the smallest eigenvalue goes to zero, 
indicating a loss of stability, precisely when we predict h — h'^. The pressure reaches a maximum at the onset of the shape 
instability. The cross (red online) indicates the pressure past h'^ is higher if the conformation were axisymmetric but this point 
is not physically realized. We have also inset graphical representations of the droplet for h/h" = 0.6,0.95,1, and 1.5. (b): 
Surface energy, E, varies monotonically with the gap height. The slope of the energy (inset) indicates that the force required 
to deform the droplet peaks aX h = h'^ past which the droplet becomes asymmetric and the amount of force required to deform 
the droplet decreases with increasing displacement. A loaded droplet would thus buckle at the transition to asymmetry. 




a) V^' 

FIG. 8. Surface Evolver simulations of the onset of the insta- 
bility in the case where the top surface is conical, x = 7r/8. 
Same notation and parameters as in Fig. 7. The smallest 
eigenvalue goes to zero, and the pressure reaches a maxi- 
mum, just past our calculated stability limit, at about h « 
1.04/i'^. Insets show graphical representations of the droplet 
for h/h" = 0.6, 0.95, 1 and 1.5. 



non- wetting scenario [23-25]. The peak load occurs at 
the critical shape, beyond which the droplet undergoes a 
limit-point buckling instability, and collapses. 

With X 7^ we are not assured that an axisymmet- 
ric instability can occur prior to contact between the 
upper and lower surface. Since we do not have a ana- 
lytical expression for the unstable point, when we must 



solve Ai?(r2,if) = given by Eq. (A16), together with 
F(A,f7,<![) = 10 and r{K,Vl,tl) = 1, with x = 7r/8 and 
9 = 37r/4, to obtain our predictions A « 1.575, VI « 1.84 
and a contact point fi w —0.93. We then compute our 
predicted critical separation h"^ = Z2 — — r2 tan x ~ 
1.18. As discussed the analytical prediction serves as 
an upper bound on the energy at the critical point, and 
we see when compared to the SE simulations, shown 
in Fig. 8, that the analysis slightly under- predicts the 
critical separation, and we obtain numerically that the 
droplet becomes asymmetric when h w 1.038h'^. 

The path from h = 1.5h'^ to h = 0.6h'^ was shown 
earlier in Fig. 5b) by dotted lines (red online), the cross 
(green online) along this path indicates where the insta- 
bility occurs in Surface Evolver. Similarly to the flat 
wall case, we find the droplet displays a drop in pres- 
sure after losing axisymmetry indicating the presence of 
a limit-point buckling instability. 



V. CONCLUSION 

In this paper, motivated by recent experiments by 
Nagy and Neitzel [23], we have used theory and com- 
putations to show that if a droplet, pinned at the bot- 
tom by a surface of finite area, if sufficiently deformed 
by a surface at the top, will always develop a shape in- 
stability at a critical compression - a result true for all 
values of the contact angle between the droplet and the 
top surface. After the critical compression, the droplet 
will then transition from a symmetric shape to an asym- 
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metric shape. The force required to deform the droplet 
peaks at the critical point then progressively decreases 
indicative of a buckling instability. If the deforming sur- 
face is flat then we predict the instability to occur when 
the apparent contact angle of the droplet at the pinned 
surface is tt, regardless of the contact angle of the up- 
per surface, similarly to past work on liquid bridges and 
sessile droplets. However, when the upper surface has 
non-trivial topology this criterion no longer holds, and 
a detailed stability analysis is carried out to predict the 
critical compression. An interesting question for future 
work would be to explore the effects of surface curvature 
on the shape instability. 
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Appendix A: Perturbations 

To satisfy the requirements on the functions F and G 
we expand them in a quarter range basis as 

F{t, <t)) = Y^ ^rn,n cos(m<^) sin (^r) , (Al) 

G(t, <t>) = Y^ cos(m(/.) sin (^t) , (A2) 

where n is odd and T = {t — ti)/At. 

Using the above definitions of F and G in Eq. (All) 

we obtain 

^(-1)^B„,„ = ^(-l)'^A„,„tanx. (A3) 

n n 

Unlike the analysis in Russo and Steen [14] we cannot use 
normal perturbations with zero magnitude at the end- 
points as we wish for the droplet to be able to have a 
variable contact point on the upper surface (fixed con- 
tact angle condition). However, as shown in Rcf. [14] the 
lowest modes are the ones with the largest increase in 
surface energy, therefore in order to establish an upper 
bound on the energy at the critical point we restrict our 
analysis to the lowest two modes in t and (j). 

We first consider the case x = 0, then generalize. With 
Eq. (A3) we have .6^,1 = -Bm,3- The boundary condition 
places no restriction on g and hence we let Am,3 = 0. 
As the surface energy is nonlinear in the shape pertur- 
bations, different modes couple and we must be careful 
in our selection of allowed perturbation shape. While 
we cannot make the perturbations normal to the shape 
everywhere, we impose a related constraint by letting 
Bm,i = ctAm,i, where a is a free parameter. This re- 
striction on the shapes reduces the unknown coefficients 
by one to facilitate the analytical calculation, while the 



free parameter allows some flexibility on the shape of the 
perturbations, and we tune its value so as to give the 
earliest possible instability. Because our choice for the 
space of allowed perturbations may not optimally mini- 
mize the surface energy, our analysis is thus only able to 
derive an upper bound for the stability limit. However, 
as we will show, this upper bound will coincide with the 
lower bound, when x = 0- With these assumptions, our 
perturbations now take the form 



F{t,cf>) = (Aoi+AnCOS0)/(i), 
G(t,</>) = (^01 +^iicos(/>)sf(i). 



where 



/(() = b,,(^t). 



g{t) = a 



+ sin I 



(A4) 
(A5) 

(A6) 
(A7) 



For a general value of x 7^ 0, we can rotate our pertur- 
bations to the X plane [/'..g']"^ = Rx[/, ,g]"^ where R,^ is 
a two-dimensional rotation operator of angle x- Hence, 
for all X) we may write 



fit) 



9it) = 



(cosx — Qfsinx) sin (2"^) ~ ccsinxsin ( ~^^ 



(A8) 



(acosx + sinx) sin + acosxsin ( 



(A9) 

We see that G2 = F2 tanx for all a. 

We require the asymmetric shapes to satisfy the con- 
tact angle condition at t2 and hence Z2 = R2 tan(^ -|- x) 
which leads directly to 



G2 



F2tan((9 + x)- 



(AID) 



Eq. (AID) merely states that our asymmetric pertur- 
bations must preserve the slope imposed by the con- 
tact angle at t2- To satisfy this for all 6 and x we set 
F2 = G2 = 0. Furthermore we require Fi = Gi = so 
that the droplet remains pinned at ti. Additionally, our 
perturbations must be directed along the surface at the 
upper bound and hence we need 



G2 = F2tanx. 



(All) 



We now expand F = ^e^F^^) and G = Y.e'G'^^^ 
where e is a small dimcnsionless parameter. We must 
ensure our perturbations, order by order, conserve vol- 
ume 

1-12 /•27r 1 ^27r 



Vl 



1 



R^Zdtd(j)- -tanx I Ridcf), (A12) 



where the second term in Eq. (A12) takes into account 
the inner cone. Expanding and subtracting off the un- 
perturbed case we get to leading order 
rt2 



Ay, = eTT^W 



/ 



r {rg + 2/i) dt - /2r2 tan x 



= 0, 
(A13) 
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hence Aj^-^ 
relationship 



0. The O(e^) term yields the foUowing 



°' ~ V 11 J 2 r (r5 + 2/i) - /ar^ tan x ' 

(A14) 

hence with A'q^ = — (^iV) where / is the quotient 

in Eq. (A14), volume is conserved to order O(e^). 
The perturbed surface energy is given by 



E =7 







7 cos 61 r^"" 
2cosx Jo 



R^d(t>. 



dtdcj) 
(A15) 



Expanding the integral and subtracting off the axisym- 
metric contribution, the leading order term of the differ- 



ence is given by 



where 



rw^ 1 / 



-=/ 

+ r' {f'+g')+{gr-fzf)^ 



1 COS0 

2 cosx 



/2(/2-47r2), 



V{t) = + 

w{t) = ff + gz. 



(A16) 



(A17) 
(A18) 



Note from Eq. (A16) that our marginal stability curve 
will be independent of 7 and A^^. We are interested in 
all outward bulging shapes, nodoids, and hence we check 
all n e [7r/2,7r). 
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